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ABSTRACT 



Aims. Our goal is to see if there is molecular gas extending throughout the optical low surface brightness disk of the galaxy Malin 2. 
Methods. We used the heterodyne receiver array (HERA) mounted on the IRAM 30m telecope to make deep observations at the 
frequency of the CO(2-l) line at nine different positions of Malin 2. With a total observing time of 1 1 hours at a velocity resolution 
of 1 1 km/s we achieve a sensitivity level of 1 mK. 

Results. We detect C0(2-l) line emission from Malin 2. The line is detected in four of the nine HERA beams; a fifth beam shows 
a marginal detection. These results not only confirm that there is molecular gas in the disk of Malin 2, but they also show that it is 
spread throughout the inner 34 kpc radius as sampled by the observations of the galaxy disk. The mean molecular gas surface density 
in the disk is 1.1 ± 0.2 Mg pc^^ and the molecular gas mass lies between the limits 4.9 x 10** to 8.3 x 10* Mg. The observed velocity 
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1. Introduction 



dispersion of the molecular gas is higher (~ 13 km s ) than in star forming galactic disks. This could explain the disk stability and its 
low star formation activity. 

Key words. Galaxies: spiral - Galaxies: kinematics and dynamics - Galaxies: ISM - ISM: molecules - Radio lines: ISM 
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Low Surface Brightness (LSB) galaxies are the most unevolved 
, class of galaxies in our nearby Universe (Impey & Bothun |1997l l. 

■ They are optically dim with diffuse stellar disks (Auld et al. 
\ 12006 ). massive HI gas disks (O'Neil et al. 2004; Matthews, van 
, Driel, Monnier-Ragaigne 2001 ) but have low star formation rates 

■ compared to regular spiral galaxies (McGaugh 1994j. They are 
halo dominated galaxies (de Blok & McGaugh 19961 Kuzio de 
Naray, McGaugh & de Blok l2008l Coccato et al. 2008 ); this 
may account for the weak spiral arms and small bar perturba- 

, tions observed in these galaxies (Mihos, de Blok & McGaugh 

■ 119971 Mayer & Wadsley I2004I I. Although the most commonly 
' observed LSB galaxies are the dwarf LSB galaxies (Sabatini et 

al. 2003 ), a significant fraction of LSB galaxies are large spi- 
rals having prominent bulges (Beijersbergen, de Blok & van der 
, Hulst il999D . These giant LSB (GLSB) galaxies have extended 
LSB disks that are poor in star formation and dust (Rahman et 
al. 120071 Hinz et al. l2007l l. The bulge dominated GLSB galaxies 
often show AGN activity (Schombert 1998, Das et al. 2009 ). 

Even though the optical properties of LSB galaxies have 
been investigated in great depth, not much is known about their 
molecular gas content. This is important as a knowledge of the 
cold gas distribution in LSB galaxies will help us understand star 
formation processes in these galaxies. Surveys of LSB galaxies 
show that they have fairly massive HI disks that may be more 
than twice the size of the optical disk (de Blok et al. 119961 
Pickering et al. 1 19971 Das et al. 2007 ). In this paper we examine 
the molecular gas distribution in a GLSB galaxy and see how it 
relates to the overall star formation in its' disk. Molecular gas 



has been detected in only a handful of LSB galaxies (O'Neil, 
Hofner & Schinnerer 2000; Matthews & G ao 120011 O'Neil, 
Schinnerer & Hofner |2P03, Matthews et al. l2005l Das et al. 
I2006 i In most cases the galaxies were large spirals with ex- 
tended optically dim disks. The low detection rate of molecular 
emission from LSB galaxies is probably due to several factors 
related to the poor star formation rate in these galaxies (e.g. de 
Blok & van der Hulst |1998l l; factors such as the lower dust con- 
tent, lower metallicity and the lower surface denity of cold, neu- 
tral gas in these galaxies. All of these properties lead to a slower 
rate of gas cooling and molecule formation. For the few galax- 
ies where molecular gas has been detected, not much is known 
about the gas extent and distribution. Such information is impor- 
tant if we want to understand star formation and disk evolution 
in LSB galaxies. 



To investigate the molecular gas and star formation in GLSB 
galaxies we studied the CO distribution in a galaxy where molec- 
ular gas has been detected, F568-6 or Malin 2 as it is widely 
known (Das et al. 120061 ). It is a nearly face-on GLSB galaxy 
at a distance of 201 Mpc. It has a prominent bulge and a very 
extended LSB disk. Its parameters are summarised in Table 1. 
There are several localized star forming regions distributed over 
its inner disk. Its metallicity is one third of the solar metallicity 
in value which is relatively high for an LSB galaxy (McGaugh 
11994 ). The CO observations of Malin 2 were conducted using 
the HERA instrument mounted on the 30 m IRAM telescope. 
Our main aim was to examine the molecular gas distribution; 
determine its extent, total gas mass and surface density. 
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2. Observations 

During March 2007 we observed the CO( 2-l) li ne in Malin 2 
with the HERA beam aiTay (Schuster et al. 120041) on the IRAM 
30m telescope at a fequency of 220.372 GHz. We specifically 
used this array as it has a wide field of view and good sensitiv- 
ity. HERA is made of 9 receivers in a 3 x 3 array spaced by 24" 
on the sky. The backend used was the Wideband Line Multiple 
Autocon-elator (WILMA). The total bandwidth was 930 Hz; it 
was made up of 465 channels of 2MHz each. The typical system 
temperatures Tsys were in the range 200-250 K for 8 receivers; 
one receiver had a systematically higher Tsys in the range 350- 
450K. The mean FWHM of each of the nine beams is 11.7". 
As our main goal was to detect the CO emission line from the 
disk of Malin 2, we kept the array fixed on the sky in the stan- 
dard pointed mode with the central beam pointed to the galactic 
nucleus; the field was tracked to get deep integrations at fixed 
positions in the disk. This enabled us to achieve good sensitivity 
but was at the cost of undersampling the source. We obtained 
deep observations at nine positions of the galaxy with a total on 
source integration time of 1 1 hours (Figure 1). The observations 
were frequency switched. The intrinsic velocity resolution was 
2.7 km s"' ; the data was then smoothed using a Hanning squared 
function. 

We reduced the data using the CLASS software of the 
GILDAS package by fitting a first order baseline to all spec- 
tra within a window going from -400 to +400 km s"' about the 
galaxys' systemic velocity of 13830 km s"'. This window was 
the same for all nine spectra. The noise level is not the same 
for all receivers; the lowest is 0.7mK and the highest is 1.6 mK 
in the 10.9 km/s channels. The noise is also not uniform across 
the band. Hence we computed the baseline and the noise in the 
-400 +400 window using emission free channels only. All the 
details are given in Table 2. The conversion factor from K to Jy 
is 8.6 Jy K ' . To convert from the antenna temperature scale 
to the main beam temperature T,„h we multiplied T* by the fac- 
tor FeffjBeff - 1 .67 where Fejf is the forward beam efficiency 
and Beff is the main beam efficiency. 

3. Results 

We have detected CO(2-l) emission from several positions in 
the disk as well as from the center of Malin 2. In the following 
paragraphs we present our results and discuss their implications. 

1. CO(2-l) Detections : Figure 2 shows the CO(2-l) emission 
spectra observed from nine locations across Malin 2. The off- 
sets from the galaxy center are indicated in each box. CO(2-l) 
emission was detected from four out of nine positions at line 
intensities above 3cr. The line at (-24, 24) is a a hint of emis- 
sion rather than a sure detection. At (0, 0) the line is broad but 
is a 3cr detection. We estimated the line parameters (flux, width 
and central velocity) by fitting a gaussian to each spectrum. We 
determined the noise level for each of the nine positions. The re- 
sults are listed in Table 2 and the best fit gaussians are overlaid 
on the spectra in Figure 2. 

2. Comparison with previous CO observations : Molecular 
gas has been detected earlier in this galaxy by Das et al. (2006) 
with the IRAM 30m telescope and the BIMA array. However, 
due to a correlator setup problem the BIMA map was incom- 
plete, it was covering a velocity range corresponding to CO(l-O) 
line emission from the east side of the galaxy only. In addition, 

' Ihttp : //www ■ iram. fr/IRAMFR/GILDAS] 



Table 1. Galaxy Parameters - Malin 2 



Parameter 


Value 


Reference 


Other Names 


F568 - 6 


NED 


Distance 


201 Mpc 


NED 


Heliocentric Velocity 


13830 km j"' 


NED 


Position (epoch 2000) 


10''39"'52'.5, +20°50'49" 


2MASS 


Size (f^j) 


1.67 


(a) 


Position angle 


70° 


NED 


Inclination 


38° 


ia) 


HI Linewidth (W^o.cojt) 


637 km i""' 


mc) 


HI Mass 


4.2 X 10'" Mq 


(b) 


Dynamical Mass 


2.5 X 10'2 Mo 


(b) 


(R < 107 Kpc) 






Disk Central Brightness 


22. 1 mag arcseC^ 


(b)(d) 


Disk Scale Length 


18.8 Kpc 


(b) 


Total R Band Magnitude 


-23.6 mag 


(b) 



" Matthews et al. (2001) 
Pickering et al. (1997) 
Corrected for inclination 
Not corrected for extinction 



the 30m telescope was pointed toward two directions only at 7" 
and ~35" east of the nucleus (see Fig. 1). Although the source 
was observed simultaneously at 1 15 GHz and 230 GHz with the 
30m telescope, only CO(l-O) emission was detected; the CO(2- 
1 ) line was not detected from either positions. Our present obser- 
vations cover the CO emission over a wider velocity and spatial 
extent than these previous observations and hence give a bet- 
ter idea of the distribution of molecular gas in the inner disk of 
Malin 2. These observations are also the first detections of the 
CO(2-l) line in Malin 2. Our detection lies below the detection 
limit of the older observations of Das et al. (2006). Their CO(2- 
1) emission spectrum had a noise (rms) of 2.6 mK whereas our 
line detections have a peak of 2 to 3 mK (Figure 2). 

Near the galaxy center, where there is overlap with the pre- 
vious CO detection, the CO line velocities are in the same direc- 
tion and are similar in shape and in width. Close to the center 
of the galaxy the CO(I-O) line as detected by Das et al. (2006) 
is broad (~200kms"') and asymmetric with the red side being 
slightly more prominent than the blue side. The CO(2-l) line 
detected at the center of the galaxy by our new HERA obser- 
vations is 243 + 76 km s"' wide and at a centroid position off- 
set by 77 ± 34 km s"' relative to the systemic velocity given in 
Table 1 . Investigating the cause of this offset would require in- 
terferometric observations. As the beams and pointing directions 
are different (23" for the CO(l-O) line and 11" for the CO(2-l) 
line) the fluxes cannot be directly compared. But if the molec- 
ular gas were uniformly distributed we would expect a line ra- 
tio in temperature scale close in the range of values found by 
Braine & Combes (1992) with the same telescope i.e., 0.4 to 
1.2. We checked this for the HERA central beam whose point- 
ing direction is the closest to the central pointing of Das et al 
(2006). A rough comparison of the brightness in mK km s^^ 
gives /co(2-i)//co(i-0) = 0.43/1.01 = 0.43 (where the CO(2-l) 
line intensity from the HERA central pointing is 0.43 and the 
CO(l-O) intensity from IRAM is 1.01). This value is at the lower 
end of the line ratio range for most galaxies but given the uncer- 
tainties (e.g. gas distribution), the value shows that the CO(2-l) 
line intensity we measure is consistent with the previous CO(l- 
0) detection. 
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3. Comparison with HI observations : High sensitivity HI 
VLA observations of Malin 2 were presented by Pickering et al. 
(1997). These reveal the distribution of the HI gas and its kine- 
matics at a spatial resolution of ~19". The HI disk is lop-sided, 
the HI emission being significantly more prominent and extend- 
ing farther out on the western side of the nucleus than on the 
eastern side. The central region is essentialy devoid of HI gas. 
Although we have only a partial view of the molecular disk due 
to our under-sampled data its properties have similitude with HI. 
Based on the HI velocity field, the centroid velocity of the HI line 
profiles at the positions (24,0), (-24,0) and (-24,24), respectively 
-h90, -150 and -110 kms"' relative to the systemic velocity, are 
in excellent agreement to those for the molecular gas (Tab. 2, col 
e). Hence the HI and molecular gas appear to share a common 
velocity assymmetry when considering the positions (24,0) and 
(-24,0). 

The higher spatial and spectral resolutions in CO allow to 
better constrain the velocity dispersion of the gaseous compo- 
nent. The CO profile with the narrowest line width is at (-24,0). 
Given the spectral response, this leads to a deconvolved velocity 
dispersion of 13 kms '. It is the best constraint that we can get 
for the velocity dispersion of the gas o-gas in the disk of Malin 
2 given the facts that a) it is obtained in a direction close to the 
major axis where beam smearing effects due to the rotation of 
the disk are minimum and b)we do not observe the blue-shifted 
wing or secondary velocity component that Pickering et al. dis- 
covered and had to blank to determine and analyze the rotation 
curve. These authors interpreted this secondary component as 
high velocity clouds, gas infalling into the disk southwest of the 
nucleus. It is best seen ~20" southwest of the nucleus where it 
is blue-shifted by as much as ~ 100 km s ' . Given the correspon- 
dance with HII regions in that direction and the fact that there 
is a relatively moderate HI column density there they suggested 
that this infall induces star formation (Wyder et al. 120091 ). We 
do not detect CO emission at (-24,-24) which is only ~6" fur- 
ther away in the southwest. It could be that the disk component 
in that region is depeleted not only in HI but also in molecular 
gas. CO observations spatially fully sampled would be required 
to confirm this hypothesis. 

In Table 2, the maximum CO(2-l) line velocities are reached 
in the offsets (24,0) and (-24,0) which is consistent with the po- 
sition angle being close to horizontal ~ 70° (2MASS images; 
Jarrett et al. 120031 . We note a velocity asymmetry between the 
emission at (24,0) and that at (-24,0), which could be due to 
infalling gas or a lopsided potential. We can estimate the rota- 
tional velocity from these lines. Assuming an inclination of 38° 
(Matthews et al. 2001| it is of the order (~ 150/ sin(0 i"' 
which implies a dynamical mass within the 24" (i.e. 23 kpc) ra- 
dius of ~ 3 X 10" Mo. 

The broad line at the center has a width close to twice the 
velocities at (24,0) and (-24,0) offsets. This could be due to the 
rotation curve turning point being inside the central beam, i.e., 
inside a radius of 6" (6 kpc). This central region which is es- 
sentially devoid of HI is where there is the strongest emission 
in CO; higher angular resolution CO observations are required 
to determine the rotation curve in that region. Would the curve 
continue to steadily decrease up to the nucleus that would point 
to a model where the gas in the central bulge region is kinemati- 
cally hot, e.g. because of gas infall to feed the AGN. 

4. Molecular Gas Mass: We derived the molecular gas surface 
densities for each line detection and the results are summarised 
in Table 2. The column densitites are derived from the line fluxes 
using N(H2) = X (1/0.9) / T „b[C0 (2 - 1)] dv; where 1(2- 
l)/I(l-0) = 0.9 (Braine & CombeslinQOll and the CO to H2 con- 



version factor is given by X = 2 x 10^° cmr^(K km i"') for cool 
virialized molecular clouds (Dickman et al. T986). The value of 
X is somewhat uncertain for LSB disks, given their low metal- 
licities and our poor knowledge of the ISM conditions in these 
galaxies. However, for a first estimate we used this value to de- 
termine the surface densities of molecular gas 1.moi in Malin 2 
(Col. f Table 2). Assuming the eight directions observed in the 
disk form a distribution of surface brightness values which is 
representative for the disk of Malin 2 up to a radius of 40", the 
mean CO(2-l) brightness is ~ 0.15 km s"'. The central re- 
gion has a beam averaged surface density that is 3 times higher 
than the mean brightness of the disk. The lower limit for M„,oi 
is ~ 4.9 X 10** Mq and is given by the sum of the four detec- 
tions in column /of Table 2 multiplied by the beam solid angle. 
To determine the upper limit for M„oh we used the noise (see 
Table 2, column a) to derive upper limits for the molecular gas 
surface densities for the regions where there were no detections 
including (-24,24) where there is a hint of emission. The upper 
limit to the total mass of molecular gas within R < 40" is 
thus M„,oi ~ 8.3 X 10^ Mq. If the molecular disk does not ex- 
tend beyond 40", then the molecular gas mass corresponds to 
~ 1.2 - 2% of the HI gas mass in the galaxy. This is similar to 
that observed in normal spiral galaxies that have molecular gas 
masses typically a few percent of the total HI mass. 
5. Star Formation threshold : The molecular gas disk in 
Malin 2 is extended and relatively massive with a significant 
amount of molecular gas but its star formation activity is low 
compared to normal spirals (Wyder et al. I2009I I. For galac- 
tic disks star formation appears to be controled by the onset 
of gravitational instabilities (Kennicutt 1989). A simple single- 
fluid Toomre disk stability model predicts threshold densities in 
agreement with the observations. This critical density is given by 
^crit = Q;/('cr^aj/3. 36 G where a: is the epicyclic frequency and 
(Tgas the velocity dispersion of the gas. From a sample of spiral 
galaxies Kennicutt found a ^ 0.69. In his study, to evaluate a, 
Kennicutt assumed that all the galaxies in his sample have ap- 
proximately the same velocity dispersion, ~ 6km i and that 
it is independent of the radial distance within a galaxy. This is a 
good approximation given what is known for the galaxies of his 
sample. Hence the prediction is that, to have disk galaxies such 
as the LSBs immune to star formation, the gas surface density 
should be below this critical density. 

Pickering et al. (1997) determined rotation curves for 4 
GLSBs. These allow them to determine the epicyclic frequency 
necessary to obtain the radial profiles of "Ecrii- They show that, by 
assuming that a - 0.69 is valid for these LSB galaxies and that 
(Tgas = 10km i"', the critical density is, in some regions, close 
or lower than the observed HI surface density, in spite of the fact 
that there is no conspicuious star formation activity there. This 
leads them to postulate that o-g^s must be larger than the assumed 
value. 

In the specific case of Malin 2 Pickering et al. find that the 
HI surface density never exceeds the critical density except at 
the location where HI is the brightest (this direction is located 
a half way between the positions (-24,24) and (-24,0)). When 
considering the 9 directions listed in table 2 only (0,0) coincides 
with an Ho- source. Using the present CO observations we can 
examine two questions: 1) does the gas surface density remain 
below the critical density when the contribution from the molec- 
ular gas component is taken into account as did Kennicutt in his 
analysis? 2) is crg„s < lOkm i"' used to determine "Lc-ii adapted 
in the case of giant LSBs such as Malin 2? Column g in Table 
2 gives the HI surface brightnesses. Comparing these with the 
molecular surface densities (col. f), the molecular gas surface 
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Fig. 1. The figure shows the R band image of the galaxy Malin 2 
with the HERA footprint overlaid as nine black circles. The 
red circles correspond to the previous IRAM observations at 
115GHz (Das et al. 2006). Each circle has a diameter equal to 
the FWHM of the corresponding beam. 



densities appear to amount to typically 1/2 to 1/3 the HI surface 
densities, with the exception for the position (-24,0) where HI 
largely dominates. Column h gives the total gas surface density 
and in col. / this quantity relative to the critical density. By us- 
ing a = 0.69, for internal consistency, the values in col. h and 
i are determined using the same hypothesis as those adopted by 
Kennicutt for the conversion factor X and the term Z which ac- 
counts for the elements heavier than H. For o-gas we adopted 
the value of 13 km s^^ that we measured with good S/N at (- 
24,0) (sect 3.3). Better constrained than the 10 km s^^ adopted 
by Pickering et al., this value is still conservative given the fact 
that it is the lowest velocity dispersion (col. d) amongst the 5 
directions where CO is detected. Formally, all the values in col. 

1 being smaller than unity it could be concluded that the disk of 
Malin 2 is sufficiently stable to prevent star formation. However, 
the directions (-24,0) and (-24,24) have surface densities close to 
the threshold preventing us to consider this as firm solid results. 
Also the observed kinematics in western and north-western parts 
of Malin 2 appears complex because of high velocity clouds 
(HVC). It could be that the star formation activity seen in Malin 

2 is more related to the interaction of these HVCs with the disk 
than to growth of density fluctuations originating from instabil- 
ities in the self-graviting rotating disk. A fully sampled image 
of the Malin 2 and high signal-to-noise to disentangle different 
velocity components and measure the line widthes are necessary 
to better understand the properties of Malin 2. 



4. Conclusions 

1. We report the first detection of CO(2-l) line emission from 
the disk of a giant LSB galaxy. Our observations reveal the 
presence of extended molecular gas in the disk of the LSB 
galaxy Malin 2. The radial extent of the molecular disk is at 
least 34 kpc and its mean line brightness ~ 0.1 5 A' km 5"' 
(on the Ti„h temperature scale). With an angular resolution 



of 11.7" this is ~ 2/3 the beam-averaged integrated line in- 
tensity in the central region. 

2. When we compared the molecular and HI gas masses of 
Malin 2, we found that the molecular gas fraction is ~1.2 - 
2% of the HI gas mass. At radii ~ 15 kpc in the directions 
where CO is detected, the molecular gas fraction is typically 
30 to 50% of the HI gas surface densities. 

3. We estimated the total surface density of the neutral gas and 
found that it is always below the critical threshold density for 
gravitational instability. Hence the disk of Malin 2 is overall 
stable against large scale star formation. 
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Fig. 2. The panel shows the CO(2-l) spectra obtained with the nine HERA receivers. The offsets in arcsecond from the galaxy center 
are marked on the top left hand corner of each figure. CO(2-l) emission was detected from five locations. The gaussian fits for these 
detections are overlaid on the spectra. The x axes are the velocity offsets from the galaxy systemic velocity (vhei = 13, 830 km s^^) 
and the y axes are the measured antenna temperatures (T^) in niK. 

Table 2. CO(2-l) fluxes, kinematical parameters, molecular, HI and total gas surface densities. 



Offset rms noise A Ico A Vfwhm <V > S„„i Eri S^as 

(arcsec) (mK) (K km s-') (km s-') (km s-') (M^pc-^) (M^pc-^) (M^pc-^) 
a b c d e f g h i 



+24 


+24 


0.7 


< 0.12 


0.39 






< 0.4 


2.1 


< 3.0 




+24 


+00 


0.9 


0.37 ± 0.07 


82 ± 20 


91 ± 


7 


1.2 ±0.2 


2.5 


4.8 


0.7 


+24 


-24 


0.8 


< 0.13 


0.45 






< 0.5 


1.8 


< 2.8 




+00 


+24 


0.9 


0.28 ± 0.05 


46 ± 10 


14 ± 


4 


0.9 ± 0.2 


2.7 


4.6 


0.7 


+00 


+00 


0.8 


0.45 ±0.11 


243 ± 76 


77 ± 


34 


1.5 ±0.4 


0.9 


3.4 




+00 


-24 


0.9 


< 0.15 


0.50 






< 0.5 


2.2 


< 3.3 




-24 


+24 


1.6 


-0.35 ±0.13 


79 ± 41 


-102 ± 


14 


~ 1.2 ±0.4 


2.2 


~ 4.4 


-0.9 


-24 


+00 


0.9 


0.20 ± 0.04 


33 ± 7 


-142 ± 


4 


0.7 ±0.1 


4.0 


5.7 


0.9 


-24 


-24 


0.7 


< 0.12 


0.39 






< 0.4 


1.8 


< 2.7 





" Aor and AS offsets relative to the position of the nucleus a, S (J2000): 10''39"'52^5, +20°50'49" 
at 10.9 km s^^ velocity resolution. 

A/co(i-o) = / Tusdv with Tmb = 1.67r^. For the non-detections, we take r^=noise and assume AV = 100 km 

fitted line widths, not deconvolved by the spectral response. 
° fitted, relative to the systemic velocity 13830 km s^' . 
*■ = 2m,,N(H2). Values for X=2.0 10^" cm'^ K km s-' (ref. sect. 3.4). 
^ from Fig. 3 of Pickering et al. (19971. 

Sgas = (1 + Z) (Shi + t^'^<^o\) cos(i). with Z=0.45, X=2.8 10-" cm-- K km i"' and i=38°. Values in brackets assume Smoi « 
' using Scrii from Fig. 10 of Pickering et al. ( il997t . re-scaled for (Tg„., = 13 fcm i"' 



